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Abstract: 

A  nematic  liquid  crystal  (NLC)  layer  with  the  anisotropy  axis  modulated  at  a  fixed  rate  q  in  the  transverse 
direction  is  considered.  If  the  layer  locally  constitutes  a  half-wave  plate,  then  the  thin -screen  approximation 
predicts  100%  -efficient  diffraction  of  normal  incident  wave.  The  possibility  of  implementing  such  a  layer 
via  anchoring  at  both  surfaces  of  a  cell  with  thickness  L  is  studied  as  a  function  of  parameter  qL  and 
threshold  values  of  this  parameter  are  found  for  a  variety  of  cases.  Distortions  of  the  structure  of  director  in 
comparison  with  the  preferable  ideal  profile  are  found  via  numerical  modeling.  Freedericksz  transition  is 
studied  for  this  configuration.  Coupled-mode  theory  is  applied  to  light  propagation  through  such  cell 
allowing  to  account  for  walk-off  effects  and  effects  of  nematic  distortion.  In  summary,  this  cell  is  suggested 
as  a  means  for  projection  display;  high  efficiency  is  predicted. 
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1.  Introduction 

Technology  of  projection  displays  exhibited  enormous  progress  in  the  past  decade.  In  particular,  Liquid 
Crystal-based  projection  Displays  (projection  LCDs),  [1,  2],  actively  compete  with  the  Micro-Electro- 
Mechanical  Systems  (MEMS)  [3],  and  with  older  technologies  like  projection  form  the  screens  of  Cathode 
Ray  Tubes  (CRTs),  [4]. 

An  interesting  subset  of  the  design  of  the  projection  LCDs  is  connected  with  the  following  idea.  In  the  “off’ 
state  of  a  pixel,  the  illuminating  beam  is  diffracted  at  a  considerable  angle.  As  a  result  the  beam  is  blocked 
by  an  aperture  positioned  downstream  of  the  light  propagation,  [5-9].  The  pixel  can  be  switched  into  the 
“on”  state  by  eliminating  the  diffraction  using  a  controlling  agent;  typically  a  locally  applied  electric  field. 
One  of  the  problems  in  the  above  design  of  the  projection  displays  is  the  requirement  of  very  strong 
suppression  of  the  intensity  left  in  the  zeroth  diffraction  order  in  the  “off”  state,  so  that  high  contrast  ratio 
may  be  achieved. 

In  this  paper  we  propose  the  use  of  an  LC  cell  with  a  specific  orientation  of  director.  Namely,  we  suggest 
using  a  cell,  with  its  director’s  azimuth  monotonously  rotating  as  a  function  of  one  of  the  transverse 
coordinates.  When  the  cell  locally  constitutes  a  “half-wave  plate”  with  respect  to  polarization,  the 
diffraction  of  incident  light  from  the  incident  diffraction  order  is  complete.  Additional  advantage  of  this 
design  of  the  element  of  a  projection  LCD  is  that  it  utilizes  both  polarizations  of  incident  light. 

In  Section  2  of  this  paper  we  describe  the  physical  principle  of  the  operation  of  the  suggested  NLC  element. 
We  develop  detailed  study  of  the  performance  of  the  suggested  display  element  in  the  rest  of  the  paper. 

In  Section  3  we  study  distortions  of  the  NLC  orientation,  i.e.  deviations  of  the  director  profile  in  the  bulk  of 
NLC  cell  from  its  ideal  structure  given  by  Eq.  1.  We  show  that  the  preferred  orientation  is  stable,  if  the  ratio 
of  thickness  to  the  modulation  period  (Li A)  is  small  enough,  typically  less  or  about  0.4.  We  also  consider 
the  change  of  this  profile  under  the  application  of  the  controlling  electric  field  (Friedericksz  transition). 
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These  results  are  obtained,  using  numeric  integration  of  the  dynamic  equations  to  find  the  equilibrium 
director  profile  of  NLC  in  the  cell. 

Then,  in  section  4  we  consider  the  problem  of  light  propagation  in  the  cell  with  the  account  of  the  walk-off 
effects.  The  latter  effects  arise  due  to  the  finite  angular  tilt  of  the  incident  and  diffracted  waves  in  the 
process  of  propagation,  so  that  the  tilted  ray  traverses  considerable  fraction  of  the  structure’s  period.  We 
estimate  the  impact  of  the  walk-off  effects  on  perfonnance  of  the  cell.  In  particular,  we  elucidate  the 
potentially  harmful  influence  of  those  effects  on  the  degree  of  suppression  of  transmitted  light  in  the  “off” 
state,  i.e.  on  the  contrast  ratio.  Finally,  we  make  a  similar  estimation  for  the  influence  of  the  bulk  distortions 
of  the  NLC  director  profile,  found  in  section  3,  on  diffraction. 


2.  Principe  of  operation  of  projection  display  element 

We  consider  nematic  liquid  crystal,  confined  inside  a  layer  0  <  z  <  L,  with  the  azimuth  of  the  director  n 
la,  b.  monotonously  rotating  in  transverse  direction  x,  Figure  1,  so  that  a  periodic  structure 

n  =  n0(jc)  =  {cos(^x),  sin(^x),  0},  q  =  2n/A,  (1) 

is  realized.  Here  A  is  the  director’s  period.  Since  NLC’s  vector  of  director  n  and  its  negative,  -n,  are 
physically  equivalent,  physical  properties  of  such  NLC  are  modulated  with  the  period  A/2.  Such  NLC  is 
locally  birefringent,  with  the  optical  axis  lying  in  (x,  v)-plane  and  rotating  along  x  direction.  Let’s  consider 
the  effect  of  this  birefringence  on  normally  incident  light  using  the  Jones  matrix  analysis. 

According  to  Jones  matrix  formulation,  the  evolution  of  electric  field  E(z)  along  the  propagation  direction  is 
given  by 


rEx{z)\  a  (Ex(0)' 
=  T(x,  z) 


(2) 


Here  E(z)  is  slowly  varying  envelope  of  the  electric  field  at  the  optical  frequency  co, 


Ereai(z,  t)  =  0.5[E(z)exp(-/cot  +  ikz)  +  compl.  conj.]. 
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Matrix  T  represents  phase  retardation  by  the  birefringent  medium  with  refractive  indexes  n0  and  ne, 
respectively.  It  should  be  written  in  local  coordinate  system  attached  to  the  director,  so  that 
A  f  cos  qx  -sin  f  exp  (i/c-z)  0  'if  cosox  sin  ax' 

T(x,z)=  .  •  .  •  ,  K  =  nine-n0)/X.  (3) 

fsingx  cosqx  J  ^  0  exp (-ik-z)j  f-singx  cosqx) 

Here  2k -z  is  the  phase  retardation  between  the  two  polarizations,  and  A  is  the  wavelength  in  vacuum.  Thus 
the  evolution  of  the  electric  field  with  distance  is  given  by 

Ex(x,  z )  =  [cos(z-k)  +  /c o s (2 t/x ) s i n (z ■  k ) ] /:A(z=() )  +  /sin(2gx)sin(z-K)iiv(z=0), 

Ey(x,  z)  =  /sin (2qx)s\ n (z ■  k )Ex(z={) )  +  [cos(z-k)  -  /co s(2t/x )s i n (z ■  k )] Z: ,(z=0 ) . 

These  relations  can  be  written  in  terms  of  circularly-polarized  components: 

R  =  EX-  iEy ,  L  =  EX  +  iEy,  (4) 

L(x,  z)=  cos(z-k)  L(z=0)  +  i'cxp(2iqx)  sin(z-K)  R(z=0) 

R(x,  z)  =  cos(z-k)  R(z=0)  +  i-cxp(-2iqx)  sin(z-K)  L{z= 0) 

If  we  assume 

L(x,  z)=  L(}(z)+  L\(z)cxp(2iqx)  and  R(x,  z)=  Ro(z)+  R  \(z)cxp(-2iqx),  (5) 

then 

Rq(z)  =  cos(z-k)  f?0(0)  +  /  sin(z-K)  ZfiO);  Zi(z)  =  cos(z-k)  Ti(0)  +  /  sin(z-K)  f?0(0).  (6) 

R-\(z)  =  cos(z-k)  R  i  (0)  +  i  sin(z-K)  To(0);  L0(z)  =  cos(z-k)  Zo(0)  +  /  sin(z-K)  R  i  (0),  (7) 

Relationships  (6,  7)  describe  the  process  of  coupling  between  two  circularly-polarized  waves,  R  and  L, 
propagating  in  the  modulated  NLC.  It  shows  that  only  the  plus  first  and  minus  first  diffraction  orders  are 
generated  by  the  device.  If  the  thickness  of  the  cell  L  satisfies  the  condition  of  half-wave  plate, 

L(ne  -  n0)  =  A/2,  i.e.  k -L  =tx/2, 

then  all  power  of  incident  wave  with  right  circular  polarization  is  transferred  into  the  (+lst)  diffraction  order, 
and  that  diffraction  order  acquires  left  circular  polarization.  Similarly,  incident  wave  with  left  circular 
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polarization  is  100%  transferred  into  the  (-1st)  diffraction  order,  and  that  diffraction  order  acquires  right 
circular  polarization, 

R-\{z=L)  =  /Lo(z=0),  Li(z=L)  =  iRo(z=0).  (8) 

This  important  property  of  the  sinusoidally  aligned  NLC  has  been  previously  discussed,  [10-12].  It 
onstitutes  the  basis  of  the  use  of  the  NLC  structure  (1)  as  an  element  of  a  projection  display.  In 
onfiguration  shown  on  Figure  2  we  propose  a  possible  design  of  such  display  cell.  Transparent  electrodes 
are  attached  to  the  cell  surfaces.  When  no  voltage  is  applied  and  the  NLC  structure  is  periodically  aligned, 
the  waves  do  diffract  and  can  be  blocked  with  an  aperture.  In  this  case  the  cell  works  in  the  “off’  regime. 
The  cell  can  be  switched  to  the  “on”  regime  by  application  of  voltage  sufficiently  high  to  destroy  the 
periodic  structure.  Then  the  NLC  cell  becomes  homogeneous,  light  not  diffracted,  and  passes  through  the 
aperture. 

In  addition  to  this  application,  this  element  can  be  used  as  an  electrically  controllable  polarization  beam 
splitter  since  the  only  diffracted  orders  are  the  two  first  orders,  which  are  right-  and  left-circularly  polarized. 
If,  however,  the  cell  is  illuminated  by  two  right-  and  left-circularly  polarized  plus-  and  minus-first  diffracted 
orders  they  combine  into  a  single  beam.  In  this  way,  this  cell  may  be  used  as  an  electrically  controllable 
incoherent  beam  combiner. 

In  the  rest  of  the  paper  we  deveploe  detailed  study  of  the  perfonnance  of  the  suggested  projection  display 
element. 

3.  Equilibrium  of  the  director  in  the  periodically  aligned  NLC  cell 

Possible  way  to  obtain  the  desired  NLC  structure  as  in  Eq.  (1)  is  through  anchoring  at  the  top  and  the 
bottom  interfaces.  Prescribed  angle  of  anchoring  of  NLC  at  the  boundary  is  usually  achieved  by  mechanical 
rubbing  of  the  surface  coated  with  a  very  thin  polymer  layer.  Orientation  of  Eq  (1)  evidently  cannot  be 
realized  with  rubbing.  However,  there  is  another  method  to  create  the  prescribed  anchoring  azimuth  which 
utilizes  photo-polymerization  under  illumination  by  a  linearly-polarized  UV  light,  [13-16].  Monotonously 
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rotating  azimuth  of  local  linear  polarization  of  the  UV  light  may  be  achieved,  if  light  constitutes  a  pair  of 
two  mutually  coherent  circularly  polarized  waves,  one  right  circular,  and  the  other  left  circular,  propagating 
at  an  angle  to  each  other.  Indeed, 

E(.x)  =  fso(ex  +  /ey)exp(-/f/x)  +  E0(ex  -  iey)cxp(iqx)  =  2ZTo[ex  cos(qx)  +  ey  sinfi/x)]  =  const-n0(x).  (9) 
The  rest  of  this  Section  3  deals  with  the  problem  of  stability  of  the  desired  structure  of  the  director,  Eq.  (1), 
in  the  bulk  of  the  NLC  cell.  When  thickness  L  is  small,  the  director  distribution  in  the  bulk  should  follow 
the  distribution  prescribed  at  the  surfaces.  However,  such  distribution  is  unstable  for  a  thick  cell,  since  its 
elastic  energy  can  be  reduced  by  formation  of  a  more  homogeneous  distribution  in  the  bulk.  The  onset  of 
this  instability  is  characterized  by  the  value  of  critical  thickness,  over  which  structure  Eq.  (1)  becomes 
impossible  to  achieve. 

In  order  to  obtain  the  director  distribution  in  the  NLC  which  is  anchored  at  the  boundaries  and  to  analyze  its 
stability,  we  have  to  solve  the  equilibrium  equations  with  appropriate  boundary  conditions.  This  solution  is 
convenient  to  achieve  via  consideration  of  time-dependent  process  of  build-up  of  the  steady-state 
distribution  of  the  director. 

Dynamics  of  the  director  of  the  NLC  can  be  described  with  the  evolution  equation,  [17,  18], 

SF  _  SR 

Sn,  k  S(dknl)  Sh ,.  ’  1 

Here  R  is  the  dissipative  function,  and  F  is  the  nematic  free  energy,  respectively: 

R  =  /  {(dn/dt)-(dn/dt)}/2,  F=  (l/2)-[Ai(div  n)2  +  A2(n-curl  n)2  +  iC(n  x  curl  n)2].  (11) 

After  substituting  (11)  into  (10)  one  can  write  the  evolution  equation  in  tenns  of  director  n 

dn/dt  =  h/y,  h  =  H-n(n-H).  (12) 

Vector  H  is  known  as  the  “molecular  field”  before  the  procedure  of  its  orthogonalization  to  local  director  n, 
while  vector  h  is  the  “molecular  field’  after  that  procedure.  Cartesian  components  of  H  are 

Hi  =  ( Ki  -  K2)di  (dana)  +  K2  ( dada )  n,  +  (K2  -  K2)  {«a»p  (<9a3p  nd  +  (dp  «p)(da  n,)  + 
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+  na  (da  «p)(3p  Hi)  -  na  (da  «p)(<3(  «p)}  (13) 

The  solution  we  are  looking  for  should  satisfy  the  boundary  conditions  at  the  surfaces  z  =  0  and  z  =  L, 

n(x,  z=0)  =  n(x,  z=L )  =  no(x)  =  {cos(^x),  sin(<7x),  0},  (14) 

In  the  general  case  this  equation,  as  written  through  the  director  and  its  spatial  derivatives,  cannot  be  solved 
analytically.  However,  an  analytical  solution  exists  in  monoconstant  case. 

Solution  and  its  stability  in  monoconstant  approximation 

In  monoconstant  approximation,  when  K\  =  K2  =  K3  =  K,  the  molecular  field  (13)  is  reduced  to 

H  =  A(V-V)n  (15) 

Direct  substitution  shows  that  the  director  distribution  n0(x)  =  {cos (qx),  sinfi/x),  0}  prescribed  by  the  Eq.  (1) 
satisfies  the  equilibrium  condition  in  the  whole  volume  of  the  cell,  since  the  H-vector  turns  out  to  be 
parallel  to  n0(x),  and  hence  h-vector  is  zero: 

3n0  idt  =  (X/y){(V-V)n0  -  n0  [n0-(V-V)n0]}  =  0, 

Therefore  no(x)  is  an  exact  solution  in  the  monoconstant  case.  The  onset  of  instability  of  this  solution  can  be 
evaluated  by  assuming 

n(x,  z,  t)  =  {cos(P)cos(gx+a),  cos(P)sin(gx+a),  sin(P)} 
where  a(x,  z,  t)  and  0(x,  z,  t)  are  small  deviation  angles.  It  can  be  approximated  with 

n(x,  z,  t)  «  n0(x)+  a(x,  z,  t)m(x)  +  ez  P(x,  z,  t),  m  (x)  =  {-sin(^x),  cos(^x),  0}  (16) 

As  a  result  one  can  obtain 

dp  idt  =  (/f/y)[c/2p  +  (V-V)p],  da/dt  =  (Ki y)(V-V)a  (17) 

It  is  convenient  to  look  for  the  solution  of  the  p-equation  in  the  form  of  a  superposition  of  eigenmodes, 
exp(-pmt)sin(nmz/L).  Then  equation  (17)  yields 

p  oc  X  |  cxp(-p,„/)sin(/«7rz/L)} ,  p,„  =  (K/y)  [(nm/L)2  -  q2].  (18) 
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Therefore  the  onset  of  the  instability  begins  (at  the  mode  m  =  1),  when  L  =  0.5A,  or  qL  =  n.  The  system  is 
always  stable  with  respect  to  the  a-perturbations,. 

Numerical  calculations. 

While  monoconstant  approximation  is  relatively  simple  and  illustrative,  it  represents  real  NLCs  only  to 
some  approximation.  To  examine  real  NLCs,  we  have  solved  the  dynamic  equations  (12,  13)  numerically. 
We  assume  arbitrary  initial  distribution  with  fixed  periodic  boundary  conditions  (14)  at  the  surfaces.  We 
then  use  Eq  (12)  to  allow  the  LC  to  relax  into  the  steady-state  equilibrium  distribution.  When  the  ideal 
periodically  aligned  structure  n0(x)  from  Eq.  (1)  is  stable,  the  volume  average 

W'fIJJ'v"'  (is) 

V  V 


is  zero.  It  usually  becomes  non-zero  in  the  presence  of  the  considerable  distortions  in  comparison  with  the 
ideal  structure.  Since  the  initial  distribution  and  boundary  conditions  do  not  depend  on  the  y-coordinate,  we 

Fig. 

3a,  b. 

Fig. 

4a  b  respectively.  It  has  a  pronounced  critical  point,  which  separates  stable  and  unstable  regions.  The  derivative 
of  (nz)  with  respect  to  parameter  (LIN)  is  presented  in  Figures  3b  and  4b;  this  derivative  has  a  sharp  peak  at 


have  omitted  y-derivatives  in  our  computation. 

Figures  3  a  and  4a  demonstrate  the  plots  of  ( nz )  as  a  function  of  thickness  for  monoconstant  FC  and  for  5CB 


Table 

1. 


the  critical  thickness.  As  expected,  critical  thickness  is  L  =  0.5A  for  monoconstant  case.  Table  1  lists  values 
of  critical  thickness  and  maximum  diffraction  angle  (calculated  using  Eq.  (26)  below)  for  a  number  of  other 
NFC’s. 

Asymptotic  analytical  solution:  thin  cell. 

When  the  thickness  is  well  below  critical,  the  sinusoidally  aligned  structure  is  reasonably  close  to  the  stable 
solution  of  the  equations  of  equilibrium.  The  director  may  slightly  deviate  fonn  the  ideal  unperturbed 
distribution  due  to  the  fact  that  the  Frank’s  constants  are  not  equal  to  each  other.  This  deviation  can  be 
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found  using  perturbation  theory  with  small  dimensionless  parameter  qL.  Assuming  n  ~  n0  +  a  m  +  ez0,  as  in 
Eq.  (16),  and  keeping  the  terms  up  to  ( qL)~  including,  one  can  get 

h0  +  K2  m  (d2a/dz2)  =  0,  h0  =  q2  (K\  -  Kf)  sin(</x)cos(</x)-m;  (3  «  0.  (20) 

That  allows  to  find  the  profile  of  azimuthal  deviation  angle  a,  if  conditions  a(at  the  boundaries)  =  0  are 
assumed: 

a(x,  z)  =  {qL)2  0.5  [(K\  -  K2)!K2]  sin(<:/x)cosfi/x)  [z(L-z)/L2]  (21) 

and  therefore 


Fig. 
5a,  b. 


6 nx(x,  z)  =  ~(qL)2  0.5  [( K \  -  K2)!K2]  sin2fi/x)cosfi/x)  [z(L-z)/L2]  (22) 

5 nv(x,  z)  =  (qL)2  0.5  \(K\  -  K2)/K2]  sin(r/x)cos2fi/x)  \z(L-z)/L2]  (23) 

Numeric  calculations  confinn  good  validity  of  this  perturbative  solution  for  small  (qL),  less  or  about  0.3. 

For  larger  values  of  that  parameter,  in  the  vicinity  of  the  critical  thickness,  actual  amplitude  of  perturbations 
was  found  to  be  much  smaller  than  the  perturbative  analytic  result.  In  particular,  for  the  NFC  5CB  and  qL  « 
2,  the  amplitude  of  nx,  //,-distortions  was  about  5  times  smaller. 

Friedericksz  transition  in  periodically  bent  NLC  under  the  applied  electric  field. 

The  pixel  transmission  is  supposed  to  be  controlled  via  application  of  external  quasi-static  (at  radio- 
frequency)  voltage  V  across  the  layer  0  <z  <  L,  i.e.  via  the  application  of  the  field  E  =  tzVIL.  Application  of 
such  field  is  accounted  for  by  the  addition  of  the  term  -0.5so(S||  -  s_i_)(E-n)  to  the  free  energy  density  F 
from  Eq.  (11).  The  analysis  analogous  to  the  one  given  above  shows  that  the  instability  threshold  is  given 
by  a  line,  which  in  monoconstant  approximation  is  an  ellipse  in  (L,  f)-planc: 

(qL)2  +  [so  (S||  -  £±)  /  K2]  V2  =  n2,  (L  /  L0)2  +  ( VI  V0)2  =  1 ,  (24) 

2  i/  —12 

where  L0=A/2  and  Vo  =  {n  K2  I  [so  (S||  -  s±)]}  2.  Here  So  is  8.85-10  “  Farad/m,  sy  and  s_l  are  the 
dimensionless  longitudinal  and  transverse  dielectric  constants,  respectively. 

Figures  5a  and  5b  show  the  result  of  numerical  build-up  of  the  steady-state;  namely  the  value  of  (nz),  as  a 
function  of  cell  thickness  L  and  applied  voltage  V  for  a  cell  with  monoconstant  NFC  and  for  5CB 
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respectively.  On  Figure  5a  the  stability  area  inside  the  ellipse  is  the  region  where  the  periodically  aligned 
structure  exists  in  ideal  form  Eq.  (1).  Our  computations  show  that  the  form  of  such  a  region  is  close  to 
ellipse,  even  if  the  NLC  in  consideration  is  not  monoconstant. 

4.  Diffraction  by  periodically  aligned  NLC 

As  it  was  shown  in  section  2,  achieving  100%  diffraction  efficiency  requires  that  NLC  layer  locally 
constitutes  X/2  plate,  i.e.  that 

L-(ne  -  n0)  =  X/2.  (25) 

On  the  other  hand,  the  diffraction  angle  can  be  found  from  the  equation  sinO  =  2X/A,  where  taking  half  of 
the  mathematical  period  A  accounts  for  the  equivalence  of  plus  and  minus  directions  of  the  director.  Thus, 

sin0  =  2k/A  =  4(ne-no)(Z/A).  (26) 

Through  this  equation  the  upper  limit  of  L/A  imposes  a  limit  on  the  maximum  diffraction  angle,  shown  in 
the  Table  1. 


Fig.  6. 


The  above  conclusion  about  100%  diffraction  efficiency  was  made  in  the  approximation  of  a  thin  phase 
screen.  In  other  words,  up  to  now  we  have  ignored  the  effects  of  walk-off  of  the  ray  across  the  periodic 


structure,  as  the  light  propagates  at  an  angle  (0mc/n)  to  the  z-direction  inside  the  cell,  see  Figure  6.  Another 
effect  which  may  deteriorate  the  performance  of  the  suggested  device  is  the  influence  of  possible  deviations 
of  the  director  from  the  ideal  structure  Eq.  (1).  To  account  for  both  of  these  effects,  we  have  to  write 
differential  equations  that  describe  1)  propagation  of  tilted  waves,  and  2)  coupling  of  different 
polarizational  /  angular  components  through  anisotropic  dielectric  pennittivity  tensor. 

Since  the  values  of  the  angle  are  relatively  small,  especially  inside  the  material  of  the  cell,  we  will  use 
paraxial  (parabolic)  approximation  for  the  dependence  of  z-component  of  wave  vector  kz  on  sin(0)  «  0,  [19]. 
We  present  the  electric  field  of  the  light  wave  inside  NLC  in  the  form 

Ereai(z,  t)  =  0.5[E(x,  z)exp(-zcot  +  ikz)  +  compl.  conj.],  k  =  (con/c),  n  =  ( ne  +  n0)/ 2,  (27) 


Efx,  z)  =  0.5(ex  +  iey)R(x,  z)  +  0.5(ex  -  iey)L(x,  z), 
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R(x,  z )  =  exp(7.s’x)  ^  Rm  ( z )  exp(2 iqmx),  L(x,  z )  =  exp(/.sx)  ^ Lm  (z)  exp(2 iqtnx),  (29) 


where  the  amplitudes  Rm(x,z)  and  Lm(x,z)  satisfy  equations 
dRm  i(s  +  2mq)2 


dz  2k 

dLm  i(s  +  2mqY 


K  0)  +  ^Rn+P(z)PP  +  iTjL>n+P  ^KP 


Lm  (Z)  +  1YjL^P  (Z)ap  +  R  >n-p  (z)(^  )  : 


(30) 


(31) 


dz  2k 

f  f 

The  variable  s  =  (co/c)0inc  characterizes  the  angle  of  incidence  of  the  input  wave  (in  air).  The  terms 
(s+2mq)  12k  in  Eqs.  (30,  31)  constitute  the  difference  between  k:  ,  z-projection  of  the  wave  vector  for  the 
given  m-th  Fourier  component  in  the  expansion  (29)  and  the  value  k  =  (co n/c),  which  was  assumed  in  Eq. 
(27)  for  normal  incidence.  To  be  more  precise,  we  took  the  paraxial  (parabolic)  approximation  for  this 
difference.  These  tenns  describe  propagation  of  tilted  beams  and,  thus,  walk-off  effects.  Coupling 
coefficients  pp,  op  and  kp  describe  the  influence  of  optical  anisotropy  of  liquid  crystal  in  periodically 
aligned  configuration. 

The  approximation  of  ideal  structure  Eq.  (1)  corresponds  to  only  one  non-zero  coupling  coefficient: 

Ki  =  (ki)*  =  k  =  n-(ne  -  noyk.  (32) 

Meanwhile  all  the  other  coupling  coefficients  pp,  <jp  and  kp  are  zero  for  propagation  in  the  ideal  structure 
Eq(l).  Then  the  system  (30,  3 1)  for  the  incident  R0  -wave  simplifies  to 

dR0(z)/dz  =  -iV()Ro(z)  +  i  k-L,(z);  dL\(z)tdz  =  i  k -Ro(z)  -  iv i  L\(z),  (33) 

Vo  =  s2/2k,  v  i  =  (s  +  2qf/2k,  (34) 

Solution  of  the  system  (30),  (3 1)  is  well-known: 

ncfj/T-t  —  /  u\  «inf  j/7-t  /  w  <i\r\(nv\  /  u  \ 

(35) 


rw 

(Rn(z  =  0)^ 

^cos(wz)  -  /(£  /  u)  sin(wz) 

i/cs'm(uz)/  u  ' 

=  M(z)  • 

o 

II 

1 1 

,  M  =exp(i(p) 

v  i/cs'm(uz)/  u 

cos(wz)  +  i(^/  u)  sin(  uz)y 

(p  =  -(v0  +  Vi)/2, 


\  =  (Vo  -  Vi)/2=  -{sq  +  q2)lk,  u  =  (k2  +  S2)0'5, 


(36) 
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Those  expressions  show  that  the  diffraction  efficiency  r|,  i.e.  intensity  coefficient  of  transformation  from  R0 
into  L], 

r)  =  [(k/w)  sin(wz)]2  (37) 

has  its  maximum  at  £,  =  0  for  any  value  of  z.  in  the  range  0  <  z  <  ti/2k.  Moreover,  100%  efficiency  is 
achieved,  when  k L  =  n-(ne  -  n0)/XL  =  %I2,  i.e.  when  the  condition  Eq.(25)  of  half-wave  plate  is  satisfied. 

In  view  of  our  task  of  creating  an  element  for  a  display  pixel,  we  are  mostly  concerned  with  leaving  no 
intensity  in  the  original,  un-diffracted  order.  This  remaining  intensity  will  constitute  the  following  fraction 
of  the  incident  intensity  \Ro\~'- 

\Ro(z  =  Lf  =  (l-ri)-|7?o(z  =  0)|2  *  | R0(z  =  0)|2-[(2^/tt)2  +(k L  -  tt/2)2].  (38) 

We  see  that  even  for  the  optimum  value  of  interaction  strength,  k L  =  n/2,  remaining  transmission  is  about 
0.4(^T)~.  In  its  turn,  the  optimum  value  c  =  0  is  achieved  at  “symmetric  geometry,  0inc  =  -(k/A).  Indeed, 

2^/n  =  -(2L/ny(sq  +  q2)/k  =  -[4  L/(nA)]-[Qmc  +  O.5|0dlft(]  =  -[4L/(nA)]-[Qmc  +  (k/A)]  (39) 

Let  us  make  an  estimation  of  the  un-diffracted  intensity  in  the  worst  case  scenario  of  the  normal  incidence 
of  illuminating  wave,  0inc  =  0: 

(l-r|)  *  [4Lk/(/;A2)]2  =  J  [8(/?c  -  n0)/n\-(L/A)2}2  (40) 

For  example,  taking  (LI A)  =  0.2  (below  critical  value)  and  ( ne  -  n0)/n  =  0.1,  we  get  (l-p)  «  1.0-1 0  3  ,  or 
l/(  1— r|)  «  1000.  This  last  number  characterizes  the  value  of  potential  contrast  ratio,  i.e.  intensity  ratio  for 
bright  state  of  pixel  to  that  for  the  dark  one.  This  estimate  shows  that  the  walk-off  effects  will  not  degrade 
the  performance  of  the  suggested  device  to  a  considerable  degree. 

Another  potential  source  of  transmission  of  incident  light  into  zero-th  diffraction  order  is  the  deviation  of 
equilibrium  structure  of  a  non-monoconstant  NLC  from  its  ideal  profile  Eq.  (1).  Taking  the  perturbative 
approximation  for  those  distortions,  Eqs.  (21-23),  one  gets  additional  coupling  coefficients  between  various 
diffraction  orders.  Suppose  for  definiteness  that  the  light,  which  illuminates  the  pixel),  is  an  incoherent 
combination  of  right-  and  left-circularly-polarized  waves.  Then  the  most  “damaging”  is  the  coefficient  kq(z), 
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which  couples  those  orders.  In  zeroth  approximation  the  amplitude  Rq(z)  =  R0(0 )  cos(kiz)  =  A0(0)  cos(kz), 
and  then  coupling  via  Ko  results  in  excitation  of  additional  amplitude  5Zo(z),  incoherent  with  the  incident 

MO); 

SL0  ( z  =  L)  =  £  k0  ( z)R0  ( z)dz  ,  (41) 

Calculations  show  that  the  distortions  given  by  Eqs.  (21-23)  lead  to  the  following  expression  for  k0(z),  and 
thus  to  the  incoherent  addition  5Co(z): 

k0(z)  =  k { [(/C  -  Ki)/  K2](2nL/A)2} [z(Z-z)/(4T2)] ;  5 L0(z=L)  =  0.043/ { [(/f3  -  A",)/  K2](2nL/A)2}R0(z=0). 

where  we  used  the  already  assumed  condition  k L  =tt/2.  Taking  the  worst  case  scenario,  when  L/A  =0.3,  i.e. 
when  thickness  is  close  to  critical,  for  5CB,  and  recalling  that  for  thickness  close  to  critical  the  numerically 
found  amplitude  of  distortion  is  about  5  times  smaller  (for  5CB)  than  the  result  of  analytic  calculation  by 
perturbation  theory  ,  Eqs.  (21-23),  we  obtain 

\8L0(z=L)\2  =  |A0(z=Z)|2-(0. 034/25)  =  |A0(z=Z)|2  /727. 

We  see  that  the  deviation  of  the  director  profile  from  the  ideal  structure  Eq.  (1)  does  not  significantly  harm 
the  contrast  ratio. 

5.  Conclusion 

To  conclude,  in  this  paper  we  suggest  a  new  design  of  an  element  of  LC  projection  display  based  on 
transverse  periodicity  of  azimuthal  alignment  of  nematic  liquid  crystal.  We  have  calculated  the  conditions 
of  stability  of  such  LC  structure  as  well  as  the  conditions  when  this  element  completely  eliminates  the 
zeroth  order  of  the  incident  light  via  diffraction.  Several  important  factors  that  might  deteriorate  the 
perfonnance  and  contrast  ratio  of  such  an  element  were  considered  in  detail.  Among  essential  advantages  of 
the  suggested  device  is  the  possibility  to  completely  utilize  the  total  power  of  depolarized  illuminating  light, 
as  well  as  rejection  of  the  signal  using  aperture  instead  of  polarizer.  Our  estimates  make  us  cautionary 
optimistic  about  the  prospects  of  creation  of  such  an  element. 
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Figure  captions  to  the  paper  Periodically  Aligned  Liquid  Crystal:  Potential  application  for  projection 
displays,  by  H.  Sarkissian,  N.  Tabirian,  B.  Park,  and  B.  Zeldovich. 

Figure  1.  Orientation  of  Nematic  Liquid  Crystal  in  a  planar  cell  0  <z  <L,  with  the  azimuth  (p(x)  of  the 
director  monotonously  growing  along  the  transverse  coordinate:  (p(x)  =  qx  =  2nxl A.  Figure  la:  top  view. 
Figure  lb:  side  view. 

Figure  2.  Design  of  a  pixel  of  a  projection  display  operating  on  periodically  aligned  NLC.  When  no  voltage 
is  applied  across  the  cell,  light  of  both  polarizations  is  100%  diffracted  into  ±l-st  orders  and  is  blocked  by 
the  aperture  of  the  projection  lens.  Application  of  voltage  results  in  suppression  of  diffraction  and  thus  in 
the  bright  state  of  the  pixel. 

Figure  3.  Instability  of  the  ideal  periodically  aligned  structure  Eq.  (1)  of  NLC  arises  at  the  thickness  L  <  aA, 
where  the  dimensionless  parameter  a  characterizes  the  threshold.  Figure  3a:  monoconstant  NLC,  for  which 
both  analytic  and  numeric  calculations  yield  a  =  0.5.  Figure  3b:  derivative  of  pervious  function  better 
illustrates  the  sharp  onset  of  instability. 
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Figure  4.  Instability  of  the  periodically  aligned  structure  Eq.  (1)  of  NLC  for  the  non-monoconstant  NLC: 
particular  values  were  taken  as  for  5CB.  Figure  4a:  symmetry-breaking  perturbation  onsets  atZ  >  0.33A. 
Figure  4b:  derivative  of  pervious  function  better  illustrates  the  sharp  onset  of  instability. 

Figure  5.  Onset  of  Friedericksz  transition  in  periodically  aligned  liquid  crystal:  volume  average  value  of 
{(nz)  )  versus  parameters  LI A  and  applied  voltage  V.  Figure  5a:  monoconstant  NLC.  Figure  5b:  graph  for  the 
5CB. 

Figure  6.  Towards  the  calculation  of  the  walk-off  effects  in  the  diffraction  process. 

Table  1.  Critical  thickness,  maximum  diffraction  angle,  and  Friedericksz  transition  voltage  for  a  number  of 
NLC’s  (input  data  was  taken  from  Ref.  [2]). 
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Figure  1 .  Orientation  of  a  Nematic  Liquid  Crystal  in  a  planar  cell  0  <z  <L,  with  the  azimuth  (p(x)  of  the 
director  monotonously  growing  along  the  transverse  coordinate:  (p(x)  =  qx.  Figure  la:  top  view.  Figure  lb: 
side  view. 
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Figure  2.  The  design  of  a  pixel  of  a  projection  display  using  periodically  aligned  NLC  cell.  When  no 
voltage  is  applied  to  the  cell,  100%  diffraction  of  light  of  both  polarizations  into  ±1 — st  orders  results  in  their 
elimination  by  an  aperture  of  the  projection  lens.  Application  of  voltage  results  in  suppression  of  diffraction 
and  thus  in  the  bright  state  of  the  pixel. 
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Figure  3  a  Figure  3b 

Figure  3.  Instability  of  the  ideal  periodically  aligned  structure  Eq.  (1)  of  NLC  arises  at  the  thickness  L  <  aA, 
where  the  dimensionless  parameter  a  characterizes  the  threshold.  Figure  3a:  monoconstant  NLC,  for  which 
both  analytic  and  numeric  calculations  yield  a  =  0.5.  Figure  3b:  derivative  of  pervious  function  better 
illustrates  the  sharp  onset  of  instability. 
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Figure  4.  Instability  of  the  periodically  aligned  structure  Eq.  (1)  ofNLC  for  the  non-monoconstant  NLC: 
particular  values  were  taken  as  for  5CB.  Figure  4a:  symmetry-breaking  perturbation  onsets  atZ  >  0.33A. 
Figure  4b:  derivative  of  pervious  function  better  illustrates  the  sharp  onset  of  instability. 
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Figure  5  a 


Figure  5b 


Figure  5.  Onset  of  Friedericksz  transition  in  periodically  aligned  liquid  crystal:  volume  average  value  of 
((nz)  )  versus  parameters  LI A  and  applied  voltage  V.  Figure  5a:  monoconstant  NLC.  Figure  5b:  graph  for  the 
5CB. 
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Figure  6. 

Figure  6.  Towards  the  calculation  of  the  walk-off  effects  in  the  diffraction  process. 
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Table  1.  Critical  thickness,  maximum  diffraction  angle,  and  Friedericksz  transition  voltage  for  a  number  of 
NLC’s  (input  data  was  taken  from  Ref.  [2]). 


LC  type 

(nm) 

An 

s// 

£1 

Ki 

( 

K2 

1012  FT 

K3 

) 

Friedericksz 

voltage 

(Volt) 

Critical 

thickness 

(A) 

Diffraction 

angle 

o 

MB  BA 

589 

0.22 

4.7 

5.4 

6.2 

3.8 

8.6 

N/A 

0.39 

20.1 

PCH-5 

589 

0.12 

17 

5 

8.5 

5.1 

16.2 

0.88 

0.33 

8.9 

PCH-5 

589 

0.11 

17 

15 

7.3 

4.5 

13.2 

2.5 

0.33 

8.3 

PCH-5 

589 

0.1 

16 

5.7 

5.9 

3.9 

9.9 

0.8 

0.35 

8 

K.15  (5CB) 

515 

0.19 

20 

6.7 

6.4 

3 

10 

0.74 

0.33 

14.7 

K21  (7CB) 

577 

0.16 

16 

6 

5.95 

4 

6.6 

0.83 

0.43 

15.6 

M15  (50CB) 

589 

0.19 

18 

6.7 

6.1 

3.74 

8.4 

0.78 

0.39 

17.5 

E7 

589 

0.23 

20 

5.1 

12 

9 

19.5 

0.96 

0.37 

19.6 

ZLI-1646 

589  0.08 

11 

4.6 

7.7 

4 

12.2 

1.2 

0.35 

6.4 
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